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Quantum cascade lasers made of nanowire superlattices are proposed. The dissipative quantum 
dynamics of their carriers is theoretically investigated using non-equilibrium Green functions. Their 
transport and gain properties are calculated for varying nanowire thickness, evidencing the transition 
from the classical to the quantum wire regime. Our calculation shows that the lateral confinement 
provided by the nanowires allows an increase of the maximum operation temperature and a strong 
reduction of the current density threshold compared to conventional terahertz quantum cascade 
lasers. 

PACS numbers: 05.60.Gg, 78.67.-n, 73.21.Cd, 73.21.Hb, 42.55.Px 



Quantum cascade lasers (QCLs) [T] are semiconduc- 
tor unipolar devices emitting coherent infrared radiation. 
Their active regions consist of stacked planar quantum 
wells (QWs) forming a superlattice along the growth di- 
rection. At first sight, a large part of the physical prop- 
erties of QCLs can be qualitatively understood in a one- 
dimensional picture, in particular for quantum confine- 
ment, tunneling processes and radiative transitions. Nev- 
ertheless, the free in-plane motion of the carriers provides 
a continuous energy reservoir (i.e. subbands) which plays 
a major role in the energy and phase relaxation processes. 
For QCLs operating in the terahertz (THz) range [2J[3], 
due to the relatively small energy separation between 
the subbands associated with the lasing levels, the non- 
radiative scattering processes are very efficient and are 
responsible for intrinsic limitations in terms of quantum 
efficiency and maximum operation temperature. 

Reducing the dimensionality of electronic systems is 
generally a way to increase the carrier lifetimes [2H7]- 
A strong magnetic field applied along the growth axis 
of a QCL quantizes the lateral motion in terms of Lan- 
dau levels. Their energy tuning with magnetic field has 
been shown to be responsible for modulation of the level 
lifetimes [51 |9J, allowing to demonstrate lasing up to 
higher temperature (225 K in ref. [TU]) than in usual 
THz QCLs (199 K in ref. pi]). However, Landau lev- 
els remain highly degenerate so that their ideal discrete 
density of states is strongly broadened by disorder effects 
|12j . On the other hand, QCLs based on 0D nanostruc- 
tures such as self-assembled or lithographically defined 
quantum dots (QDs) have been proposed [13H16I in or- 
der to achieve low threshold current and high operation 
temperature. However, a clear comparison of the device 
performances with respect to conventional THz QCLs is 
still lacking, and the intermediate regime between cou- 
pled 0D QDs and coupled 2D QWs remains to be inves- 
tigated. 

In this letter we propose a novel type of QCLs 
based on nanowire (NW) superlattice heterostructures. 
Within the formalism of non-equilibrium Green functions 



(NEGF), we present a theory of dissipative quantum dy- 
namics in wire axial superlattices, which provides a uni- 
fied model of electron transport and THz gain in super- 
lattices with and without lateral dimensional quantum 
effects. We investigate the modifications in the scatter- 
ing mechanisms and in the nature of charge transport 
associated with the transition from the 2D to the 0D lim- 
its for the lateral motion. Our calculations demonstrate 
that below a critical thickness NW THz QCLs have the 
ability to achieve higher temperature operation and lower 
current threshold than conventional THz QCLs. 

The active region of the proposed QCL consists of an 
array of parallel NWs, each of them containing a super- 
lattice heterostructure along the growth axis. The ex- 
perimental realization could be done cither by etching of 
a standard QCL heterostructure or by a bottom-up ap- 
proach. The electronic structure is modeled within the ef- 
fective mass approximation. The NW surface is assumed 
to be perfectly passivated so that no charge depletion 
or accumulation effects occur. We consider an electron 
wave function basis set of the form >f an — ip a (z)<p n (p, 9), 
where </>„ are the solutions of the Schrodinger equation 
on a homogeneous disc (the assumed NW cross-section) 
with eigenvalues E n and the ip a form a periodic basis set 
of the low energy states along the z superlattice axis [T7] . 

The electron dynamics is calculated within the NEGF 
formalism in which the perturbative scattering processes 
arc treated within the self-consistent Born approxima- 
tion (SCBA). Self-consistent (SC) NEGF implementa- 
tion of dissipative quantum transport in semiconduct- 
ing heterostructures [T8H2T] allows a coherent and non- 
Markovian treatment of the charge dynamics, where the 
scattering self-energies can be conveniently treated on a 
microscopic basis. In the present work, a SC-NEGF for- 
malism is further required since dimensionality transition 
is intrinsically a self-consistent problem between scat- 
tering processes and line broadening effects. The elec- 
tron lesser Green function (GF) reads G< p n (ti,t 2 ) = 

( a /3 n(^)a a ,n(ti)) , where a + (a) are the electron creation 
(annihilation) operators, a and f3 refering to the ip(z) 
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modes and n to the 4> n (r,9) disk eigenstates. Thus co- 
herences (i.e. a / (3) are considered in the growth direc- 
tion while only populations are considered for the lateral 
modes. Field-periodic boundary conditions [19] are ap- 
plied. 

We take into account the scattering mechanisms that 
are known to be important in QW QCLs heterostruc- 
tures: longitudinal optical (LO) and acoustical phonon 
emission/absorption, interface roughness, charged impu- 
rities and alloy disorder. In addition, scattering mecha- 
nisms specific to NW and/or quasi OD systems are also 
considered. The NW surface roughness is taken into ac- 
count as an additional elastic scattering mechanism. NW 
phonon confinement effects are considered and we calcu- 
late the coupling of electrons to longitudinal optical (LO) 
and surface optical (SO) phonons 22j, as well as confined 
acoustic phonons. 

In quasi OD systems, early studies have predicted a 
phonon bottleneck as soon as discrete electronic transi- 
tions are detuned from the optical or acoustical phonon 
frequencies |H [55] ■ More recent studies have demon- 
strated that (i) direct LO-phonon emission/absorption 
processes are totally suppressed even in case of electron- 
LO-phonon resonances due to the strong coupling regime 
between electrons and LO-phonons [24] , and that (ii) the 
dominant energy relaxation of terahertz excitations is the 
polaron anharmonicity, i.e. an indirect scattering pro- 
cess involving electron-LO-phonon polar interaction as 
well as anharmonic couplings of LO-phonons with other 
phonon modes [Jj [25] . In the present case, we have 
checked that the SCBA very well reproduces the polaron 
formation in the limits of small QDs superlattices (i.e. 
for the thinner diameter considered below), as reported 
in ref. [3U] for self-assembled QD superlattices. The elec- 
tronic lesser self-energy due to phonon interaction reads 



{t) = Y J u m D<{t)G<{t)u; 



(1) 



where U m is the interacting potential induced by the 
phonon mode m (LO or SO) and G < the electron lesser 
GF. The optical phonon anharmonicity is included in 
Eq. (JlJ by considering an interacting phonon Green func- 



tion D< instead of a free one L>„ (0) as usually done shown on Fig. |l(a)| Fig. |l(b)| shows the current density 
[HI EH HI] • However, beyond a phenomenological broad- 
ening of the optical phonon GFs [50], we calculate the 
phononic self-energies of the optical phonons due to their 
anharmonic coupling with the other phonon modes 



The elastic scattering mechanisms induced by dis- 
order effects (here interface roughness, surface rough- 
ness, alloy disorder and charged impurities) are usually 
treated within the SCBA in planar QW heterostruc- 
tures PH HH] [5T]. However, the SCBA only accounts 
for an incoherent treatment of the disorder effects. Here, 
as the NW gets thinner, the elastic scattering becomes 
mainly coherent, as the disorder potential fluctuations 
become larger than the phonon-limited broadening. In 
the present work, a coherent treatment of the disorder po- 
tential V within each lateral channel V nn (z) = (4>n\V\(f> n ) 
is achieved. For this purpose, we generate random 
correlated potentials V nn (z) from the principal compo- 
nent analysis of the covariance matrix (V nn (z)V nn (z')) . 
On the other hand, the elastic scattering between dif- 
ferent lateral eigenmodes remains treated within the 
SCBA and corresponding lesser/retarded self-energies 
read £<•*(£) = £ m#n <|K m (z)| 2 )G<' fi (£). The trans- 
port and optical properties are calculated for each ran- 
dom potential realization and finally averaged. 

We consider here the Ino.53Gao.47As/GaAso.51Sbo.49 
material system lattice matched to InP, in which a con- 
ventional THz QCL operating up to 142 K has recently 
been demonstrated 2(i . This system could be used to 
realize NW QCLs in a top-down approach, and bottom- 
up approaches may also be possible with similar mate- 
rial systems [571 HE] • We consider below 3- well resonant 
phonon designs (A and B), similar to the one holding op- 
eration temperature records in conventional THz QCLs 
[TTj . The layer sequence in nanometers of the design A is 
14/1.5/11.2/3.5/22/3.5 with barriers indicated in bold 
fonts and the doped layer underlined. The equivalent 
areal doping density per period is set to 3 x 10 10 cm' 2 . 
The interface and surface roughnesses are both assumed 
to have a rms fluctuating amplitude of 3 A with correla- 
tion lengths of 8 nm. The filling factor of the NW array 
is assumed to be unity for simplicity in the treatment 
of Coulombic interactions and for clear comparison with 
the conventional QCL geometry. 

Fig. [T] shows the calculated transport properties of the 
NW QCLs. The current-voltage characteristics calcu- 
lated for various NW diameters at room temperature are 



as a function of the NW diameter for a fixed applied bias 
of 65 mV per period for various temperatures. For suffi- 
ciently thick wires the calculated current densities satu- 



ia<(t) = J2 \W m , a , \ 2 D<°\t)D<°\t), (2) 

a,p 

in which the cubic interaction W m , a ,p follows from the 
same microscopic modeling of crystal anharmonicity as 
reported in Ref. [25j . The interacting optical phonon 
GFs .D< are then computed using the Dyson (D^ = 
- DT^ 1 - and Keldysh (D< = D^a<D^) 
relations. 



rate on Fig. 1(b) and all the IV curves for diameters larger 
than ~150 nm (not shown on Fig. 1(a) ) overlap. This cor- 
responds to the 2D regime of the lateral motion, i.e. when 
the lateral quantization effects are negligible. When the 
NW diameter is reduced below a temperature-dependent 
critical diameter, variations of the current density are ob- 



served on Fig. 1(b) corresponding to the appearance of 
lateral dimensional effects. In this quantum wire regime, 
the I-V characteristic possesses more distinct peaks as 
the reduction of the scattering processes allows adiabatic 
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FIG. 1: (Color online). Transport properties of the NW 
QCLs. (a) Current density as a function of the voltage drop 
per period for various NW diameters at lattice temperature of 
300 K, for the design A (solid and dotted lines for respectively 
stable and unstable operating regions) and for the designs B 
(dashed-dotted lines). The designed operating bias (65 mV) 
is indicated by a vertical dotted line. Layer sequences of de- 
signs A and B are explicited in the text, (b) Current density 
as a function of the NW diameter for fixed operating bias 
voltage at lattice temperatures of 100 K (red lines) and 300 
K (black lines) for the design A (solid lines) and the designs 
B (dashed-dotted lines). The top axis indicates the energy 
separation between the NW lateral ground state 4>o and the 
first excited state d>i . 



transport over larger distances. In fact, in the 2D regime 
the lateral motion provides an energy reservoir towards 
which the excess axial energy is easily transfered via the 
disorder potentials. As the lateral motion gets quantized, 
this elastic energy exchange tends to be suppressed. As 
the remaining scattering processes are mainly due to op- 
tical phonons, a strong current peak remains for potential 
drop per period matching the optical phonon energy (34 
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FIG. 2: (Color online). Intensity gain spectra at 300 K calcu- 
lated for various NW diameter for the design A (solid lines) 
and the design B300 (dashed-dotted lines). 



meV), whereas the current is strongly reduced for other 
bias voltage. For thin diameters, as this parasitic current 
peak at 34 mV becomes larger than the current at the 
laser operating bias (around 65 meV), electrical instabil- 
ity can occur in the design A. In order to restore electrical 
stability, designs B differing from the design A mainly by 
twice thicker injector and extractor tunneling barriers (7 
nm instead of 3.5 nm) are investigated. Owing to the 
sharper tunneling resonances and the effective mass tem- 
perature dependence, two slighty different designs B300 
(B100) with layer sequence in nm 14.4/1.5/11.5/7/23/7 
(14.5/1.5/11.6/7/22.5/7) are s tudie d at respectively 300 
K and 100 K. As shown on Fig. 1(a) at 300 K, the weaker 
tunneling couplings in designs B allow a reduction of the 
current leakage peak at the optical phonon energy. A fur- 
ther insight into the nature of transport is provided by 
the analysis of the spectral function (not shown here). 
For the thick injector and collector barriers in the design 
A, anticrossings in the spectral function at the tunneling 
resonances are observed only for NW diameters smaller 
than ~ 35 nm (~ 60 nm) at 300 K (100K respectively), 
indicating than the nature of these resonant tunneling 
processes evolves from an incoherent to a coherent regime 
as the NW thickness is reduced. 

We calculate the THz gain in linear response with a 
self-consistent treatment of the induced self-energies (i.e. 
including vertex corrections) [29] , which has been shown 
to be necessary to correctly treat the correlation in the 
scattering potentials [3D]. The calculated intensity gain 
spectra at the optimal operating bias are shown on fig. [2] 
at 300 K, and the peak values of the gain are plotted 
on Fig. [3] as a function of the NW diameter at 100 K 
and 300 K. The gain of the design A in the QW limit 
reaches 6 cm -1 at room temperature, which is below the 
typical reported cavity losses |31j . The gain strongly in- 



4 



E,-E (meV) 

100 10 1 



400 
200 
V E 100 









ft*'- 

t 








100 K 






/IF 




300 K 


i 





1 , 1 — I — I — I — I 1 — 

40 60 100 200 

Nanowire diameter (nm) 



on the THz QCL device performances. We have shown 
that NW confinement in the active region of THz QCLs 
is a promising way to achieve room temperature opera- 
tion. In addition, the increase of the coherent propaga- 
tion lengths in the quantum wire regime opens the route 
to designs with larger tunneling barriers and much lower 
current threshold. 
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Humboldt foundation and the Austrian Science Fund 
FWF (SFB IR-ON). S. Birner, M. Branstetter, C. 
Deutsch, P. Greek, G. Koblmiiller, M. Krall, T. Kubis, 
S. Rotter, K. Unterrainer and P. Vogl are gratefully ac- 
knowledged for fruitful discussions. 



FIG. 3: (Color online). Maximum intensity gain as a function 
of the NW diameter at 300K (black lines) and 100K (red lines) 
for the design A (solid and dotted lines for respectively elec- 
trical stable and unstable gain regions) and for the designs B 
(dashed-dotted lines, only for NW diameters satisfying elec- 
trically stability). 



creases as the NW diameter decreases into the quantum 
wire regime. Assuming cavity losses equivalent to an ab- 
sorption of 15 cm -1 , our calculations predict THz lasing 
at 300 K for NW diameters smaller than 70 nm. As the 
NW diameter is further reduced below 38 nm at 300 K 
(62 nm at 100 K), the design A becomes electrically un- 
stable (dotted lines on Fig. [3]). The design B is stable 
down to thinner NW diameters and its gain for NW di- 
ameter between 30 and 36 nm reaches values around 60 
cm -1 at room temperature, which is equivalent to the 
calculated value at 100 K in the 2D limit. Below 20 nm 
diameters, as the lateral extent of the electron wavefunc- 
tions becomes comparable with the correlation length of 
the interface and surface roughnesses, the gain strongly 
decreases due to the large energy fluctuations and the 
consequent misalignment of the QCL levels. Interestingly 
the calculated gain for the designs B300 and B100 is pos- 
itive only below a critical diameter (38 nm at 300 K and 
100 nm at 100 K), and is negative beyond up to the 2D 
limit. Indeed, coherent propagation through their large 
injector and collector tunneling barriers occurs only in 
the quantum wire regime, whereas in the 2D limit the 
tunneling is mainly incoherent and not efficient enough 
to allow population inversion. A further advantage of the 
designs B300 and B100 in the quantum wire regime resides 
in their much smaller current densities threshold (more 
than one order of magnitude smaller) than the design A 
operating in the 2D limit. 

In summary, we have developed a microscopic dissipa- 
tive quantum transport theory of electrons in NW su- 
perlattices. We have investigated the effects of lateral 
dimensionality on the physics of electron dynamics and 
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